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Pseudobinary heterostructural alloys of ZnO with MgO or CdO are studied by composing the system lo- 
cally of clusters with varying ratio of cations. We investigate fourfold (wurtzite structure) and sixfold (rocksalt 
structure) coordination of the atoms. By means of density functional theory we study a total number of 256 
16-atom clusters divided into 22 classes for the wurtzite structure and 16 classes for the rocksalt structure for 
each of the alloy systems. The fraction with which each cluster contributes to the alloy is determined for a given 
temperature T and composition x within (i) the generalized quasi-chemical approximation, (ii) the model of a 
strict-regular solution, and (iii) the model of microscopic decomposition. From the cluster fractions we derive 
conclusions about the miscibility and the critical compositions at which the average crystal structure changes. 
Thermodynamic properties such as the mixing free energy and the mixing entropy are investigated for the three 
different statistics. We discuss the consequences of the two different local lattice structures for characteristic 
atomic distances, cohesive energies, and the alloys' elasticities. The differences in the properties of Mgj:Zni_ vO 
and CdtZni_fO are explained and discussed. 
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PACS numbers: 61.66.Dk, 64.75.+g, 71.15.NC 

I. INTRODUCTION 

Recently, ZnO has attracted renewed attention as a closely 
lattice-matched substrate for GaN and also as a potentially 
useful active optoelectronic material in its own right When 
doped with Al it plays a role as a transparent conducting oxide 
in solar cells^, but also microelectronic and sensory devices 
may be fabricated involving oxides like ZnO^'^. Further, it 
is environmentally friendly, biocompatible, and tends to form 
nanostructures^. 

The tailoring of material properties is an important reason 
for combining ZnO with other group-II oxides in alloys and 
heterostructures. Mixed Mg^Zni_;rO crystals allow to tune 
the energy gap from about 3.4 eV of wurtzite {wz) ZnO to- 
wards the ultraviolet spectral region with at least 4.4 eV in 
MgvZni_vO^"''. Conversely, pseudobinary CdvZni_vO al- 
loys tend to close the gap and are, therefore, suitable candi- 
dates for optoelectronic devices in the visible spectral ranged. 
Moreover, quantum-well and other heterostructures based on 
the combinations Mg;rZni_^0/ZnO and ZnO/Cd;rZni_jO al- 
low an additional tailoring of electronic and optical properties 
by means of quantum-confinement effect a^'^^'^^ . 

The combination of these group-II oxides in alloys and het- 
erostructures raises several fundamental questions, since there 
is a crystal structure mismatch: While some binary II- VI com- 
pounds occur in the fourfold-coordinated wz structure (e.g. 
ZnO), others show the sixfold-coordinated rocksalt (rs) struc- 
ture (e.g. MgO, CdO) under ambient conditions'^ '^. Isova- 
lent and isostructural alloys of II- VI constituents are generally 
thermodynamically unstable, because their mixing enthalpy 
in either the rs structure or the wz structure is positive'^. A 
thermodynamical miscibility only exists at very high temper- 
atures, i.e., where the mixing entropy is sufficiently large; at 
lower temperatures there is a tendency for a phase separation 
of the alloys'^. However, isovalent but heterostructural II- VI 
alloys seem to be stable under certain conditions: XxZni_;tO 



alloy films occur in wz structure for x < 0.55 (X = Mg)i^"— 
or X < 0.32 (X = Cd, grown under non-equilibrium condi- 
Depending on the composition x the film prepa- 
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ration may give rise to a non-uniformity of pseudobinary thin 
films as observed experimentally for MgxZni_vO^'. With 
increasing molar fractions x (Mg or Cd content) a change 
of the coordination (fourfold to sixfold) of the atoms is ex- 
pected. Both hexagonal and cubic MgvZni_,0 thin films 
are reported^ The situation in pseudobinary CdvZni_,0 
is less clear. The rs structure of CdO seems to limit the 
equilibrium solubility in wz-ZnO to molar fractions below 
X = 0.0 7^^'^^ . On the other hand, metal-organic chemical va- 
por deposition (MOCVD) techniques produced alloys up to 
X — 0.7'^. Molecular beam epitaxy (MBE) leads to epilayers 
with a structure close to the wz one without an indication for 
a phase separation up to a Cd concentration of x = 0.32'^. 

Summarizing, there is no microscopic picture of how the 
heterostructural aspect influences the properties of isovalent 
pseudobinary alloys with oxygen as anion. Stability, sol- 
ubility, and tendency of decomposition of systems such as 
MgvZni_;rO or CdvZni_;rO are barely understood, so is their 
different behavior One reason is that most of the existing 
theoretical studies are limited to ordered structures such as 
MgZnO i^^'^^ or somewhat more complex geometries'^. 

On the other hand, the theoretical methods for the de- 
scription of the thermodynamic and structural properties of 
pseudobinary alloys of the type AvBi_^C have been devel- 
oped over the last 20 years. A basic approach is to simu- 
late the configurational average by representing all possible 
occurring local bonding configurations using certain crystal 
geometries. There is one study^^ in which the chemically dis- 
ordered MgiZni_xO alloys have been simulated using spe- 
cial quasi-random structures'^. In another approach'^'^" the 
mixing free energy of a random alloy is described combining 
a cluster-expansion method^ '"^^ with the generalized quasi- 
chemical approximation (GQCA)^^. In several publications 
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about pseudobinary alloys with the first-row elements N or O 
as anions^^'^"'^^, this method has been applied together with 
an ab-initio description of total energies and atomic geome- 
tries. However, in all these studies the (local) crystal structure 
has been fixed being wz, zb, or rs. Only recently the isova- 
lent but heterostructural MgO-ZnO system has been studied 
by combining the rs with the wz crystal structure within a 
cluster expansion^^'^^. Energetic stability has been found in 
the sixfold-coordinated structure for Zn concentrations below 
67%, giving rise to spontaneously ordered alloys"'. Without 
a configurational average and studying only certain clusters 
with fixed molar fractions^^ no conclusions are possible for 
different growth situations. Coherent potential approximation 
calculations^^ for various crystal phases allowed for the con- 
struction of a phase diagram with a phase transition from wz to 
rs close to x = 0.33. However, the coherent potential approxi- 
mation, as a single-site theory, is expected to be more reliable 
close to the binary end components, whereas for intermediate 
compositions x the single-site approximation is questionable. 

In this paper we extend the generalized quasi-chemical ap- 
proach for one local crystal structure to alloys that consist of 
two binary compounds with two different equilibrium crystal 
structures, i.e., rs and wz- The results of the GQCA are com- 
pared to two other cluster statistics, that may be considered as 
limiting cases of the GQCA with respect to the number of pos- 
sible atom arrangements. The theory is applied to heterostruc- 
tural Mg;tZni_jO and CdvZni_iO alloys. In Section HIl we 
present a cluster expansion for the rs crystal structure along 
with the combined statistics that deals with wz and rs clusters 
simultaneously. The computational approaches that are used 
for carrying out the actual calculations are explained in Sec- 
tion Hn] Results for the thermodynamic properties are given 
in Section |IV] and for the structural and elastic properties in 
Section[Vl Finally, Section|VT]concludes the paper. 



II. ALLOY STATISTICS AND THERMODYNAMICS 
A. Cluster expansion for wurtzite and rocksalt structures 

We study pseudobinary alloys A^Bi_,C with atoms 
of type C on the anion sublattice and atoms of type A 
or B on the cation sublattice. Within a cluster expansion 
method22.-i^ the macroscopic alloy is divided into an en- 
semble of M clusters consisting of 2n atoms (n anions and n 
cations) each. The total number of cations or anions is then 
given by = nM. 

From combinatorics it follows that, for a given crystal struc- 
ture, there are 2" different possibilities of arranging A- or B- 
type atoms on the n cation sites of one cluster (the occupation 
of the anion sublattice is fixed). Due to the symmetry of the 
crystal lattice, the clusters can be grouped in / + 1 different 
classes, with J depending on the actual crystal structure. Each 
class j {j = Q,...,J) contains gj clusters of the same total en- 
ergy Ej, with the degeneracy factors gj fulfilling the relation 

To each macroscopic alloy one can assign a cluster set 
{Mo, Ml ,Mj} which describes how many clusters of each 



class occur in the alloy. A single class j contributes to the 
macroscopic alloy with its cluster fraction xj that is defined 
by Xj = Mj/M. The xj fulfill the constraint 

L-^. = i, (1) 

which stems directly from the relation M — Y.j for the clus- 
ter set. The n cation sites of each cluster are occupied with nj 
atoms of species A and {n — nj) atoms of species B. Since 
the molar fraction of A atoms for the entire alloy AvBi_j^C 
is fixed by x, the cluster fractions xj have to obey the second 
constraint 

J 

^ rijXj = nx. (2) 

;=o 

Using such a cluster expansion any macroscopic alloy can be 
built from the microscopic clusters, each of which contributes 
with its cluster fraction. Consequently, within this frame- 
work each property P of the macroscopic system can be traced 
back to the respective properties Pj of the clusters. Given the 
weights Xj{x, T) for an alloy of a certain composition x at a 
temperature T and the values Pj of the property for each clus- 
ter, one can calculate the property P{x, T) for the alloy using 
the Connolly-Williams method^2i2^, 

J 

P{x,T)='£xj{x,T)Pj. (3) 

j=o 

With this approach structural, elastic, and thermodynamic 
properties, including lattice parameters and bulk moduli can 
be accessed. 

Both the Mg^Zni_;fO alloy and the Cd^Zni_^0 alloy con- 
sist of two monoxides, which show different equilibrium lat- 
tice structures: MgO and CdO crystallize in the cubic rs 
phase, while the equilibrium structure of ZnO is the hexag- 
onal wz phase under ambient conditions. In this work we 
follow previous argumentations^^'^^ that it is sufficient to in- 
clude next-nearest-neighbor correlations to capture large parts 
of the physics involved in the problem. Consequently, we use 
16-atom cluster cells (i.e., n = 8) and assume that possible 
correlations between the different clusters are small. 

While we rely on the cluster expansion for the wz structure 
which is described in Ref. 35, we present a new one for alloys 
with rs crystal structure in this work. The corresponding unit 
cells and the labeling of the cation positions for both types, 
wz and rs clusters, are illustrated in Fig. [T] In Table U we 
give the number rtj of A cations, the degeneracy factors gj, 
and one representative of the symmetry-equivalent clusters for 
each class j (characterized by the cation sites that are occupied 
with A atoms). The degeneracy gj for each class as well as 
the total number of classes 7 + 1 do not only depend on the 
number of atoms in the clusters, but also on the actual point 
group of the crystal. For n = 8 we obtain 22 classes in the wz 
case and 16 classes for the rs structure (see Table HI. Taking 
into account the degeneracies gj a total of 256 clusters per 
crystal structure are studied for each pseudobinary material 
system. 
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FIG. 1. (Color online) Illustration of atomic sites in the 16-atom 
clusters of wurtzite (a) and rocksalt (b) structure. Anions (C atoms) 
are depicted in blue (small), cations (A or B atoms) in red (large). 
White atoms represent the surrounding lattice and are not part of the 
16-atom clusters. The denotation 1 - 8 of the cations is used in Ta- 
ble|I] 



TABLE I. Cluster classes for 16-atom cells in wz and rs crystal struc- 
ture. For each class j the number jij of A cations, the degeneracy gj 
of the class, and the cation sites occupied with A-type atoms (for one 
representative of the class) are given. 
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B. Generalized quasi-chemical approximation 

In the framework of the generalized quasi-chemical approx- 
imation (GQCA),— the cluster fractions are determined 
by a minimization of the Helmholtz free energy F{x,T). We 
restrict the discussion to the mixing contribution AF to the 
free energy 



AF{x, T) = AU{x, T) - TAS{x, T). 



(4) 



The mixing contribution to the internal energy AU {x, T) is 
calculated as the sum over the contributions from the M clus- 
ters and is referenced to the value of U of an alloy consisting 
only of the two binary end components AC and BC, 



AU{x,T) 




[xej + {l -x)eo] 



(5) 



Here we introduce the definition of the excess energy Ae^ for 
the class j. 



Ae, = e^- -^ej- 



-eo 



(6) 



For the calculation of the free energy, Eq. (|4|, or thermody- 
namic properties derived thereof, an expression for the config- 
urational (or mixing) entropy has to be found. To evaluate 
the Boltzmann definition of the entropy, AS{x,T) = k^lnW, 
one has to give an expression for the number of possible con- 
figurations W. Given a cluster expansion and, therefore, a 
set of cluster fractions {xj } that fulfills the constraints ([T]) and 
(|2]l, W describes the number of possible atomic configurations 
in the entire alloy, i.e., W counts all possible ways of arrang- 
ing the A^A A atoms and A^b B atoms for one given set {jc, } 
on the ~ Na +Nb cation sites. To determine W, the num- 
ber of ways of arranging the Mq,Mi , . . . ,Mj clusters to form 
the alloy, Ml/Y[jMjl, has to be multiplied by the number of 
possibilities to arrange the cations in each cluster. Since one 
cluster of class j can be occupied by cations in gj ways, all 

M 

Mj clusters lead to gj ' possibilities. Taking into account all 
classes j one ultimately obtains 



W = 



Ml 



Us/ 



(7) 



Using this expression for W and the definition xj = Mj/M 
one finds for the mixing entropy AS in the Stirling limit 



AS(x,T) = -kBMy Xj\n \ 



(8) 



Introducing the ideal cluster fractions of a strict-regular so- 
lution^l (see Section HTcT i. 



x^ = gjx^y(i-x)"-''r 



(9) 
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the mixing entropy can be rewritten int( 



,30.34 



AS{x,T) = -kB\^N[xlnx+{l -x)ln(l ~x)] 
+ Mtx,lnQ)}. 



(10) 



Ahhough the expressions for W according to Eq. (|7]) and the 
one derived in Ref. 30 differ, Eq. ( fTOl i shows that they both 
lead to the same entropy in the Stirling limit. 

Equations (|5]l and (O fully determine the Helmholtz mix- 
ing free energy as a function of x and T, given that the 
cluster fractions xj are known for x and T. In the GQCA 
the Xj are determined by the requirement that AF{x,T) as- 
sumes a minimum with respect to the cluster distribution, i.e., 
dAF {x,T)/ dxj = 0. Hence, the Lagrange formalism with the 
constraints ([T]i and ^ yields 



GQCA 



-/5A£y 



(11) 



whereas j3 = l/k^T. The parameter Tj has to be determined 
from the condition that the x^*^*^^ obey the constraint 

The cluster fractions jc^*^*^^ according to Eq. ( fTTT l describe 
the probability for the occurrence of each cluster class /' in 
an alloy which has been prepared under thermodynamic equi- 
librium conditions that minimize the free energy. 

Since we study two cases, isostructural {wz or rs only) as 
well as heterostructural (combined wz and rs) alloys, we dis- 
tinguish three situations: (i) wz clusters only, (ii) rs clusters 
only, or (iii) simultaneous occurrence of wz and rs clusters. 
Depending on the clusters included in the calculation, the 
value for J has to be set to 7 = 21, 15, or 37 (see TablelJi and 
the minimization of AF has to be performed for given x and 
T accordingly, which leads, via Eq. (fTTI) . to the corresponding 
Xj for each case. 



C. Strict-regular solution and microscopic decomposition limit 

In the literature several techniques for the deposition of 
Mg^Zni_vO films are reported, e.g. pulsed laser deposition 
(PLD) with growth temperatures of 950... 1050 KJ°, radio 
frequency magnetron sputtering at 700 K^, and reactive 
electron beam evaporation (REBE) at a substrate tempera- 
ture of 550 K^^ Similar techniques have been applied for 
Cd^Zni_^0 layers, however, at much lower substrate tem- 
peratures, e.g. MBE with a growth temperature as low as 
450 K^^, (plasma-enhanced) MOCVD at a growth tempera- 
ture of 625 Ki^, or PLD at 700 K^. Subsequent repeated 
temper steps are reported in some cases. The films are typi- 
cally deposited on c-plane (0001) sapphire, but also a-plane 
sapphire or more exotic substrates have been used. Therefore, 
the respective experimental situation suggests the studying of 
certain non-equilibrium preparation conditions, for which the 
actual cluster statistics may be modified by kinetic barriers, 
frozen high-temperature states, as well as interface or surface 



influences. In order to simulate a dependence of the cluster 
distribution on the preparation conditions we study two limit- 
ing cases: 

(i) The strict-regular solution (SRS) modeP"*. In this case, 
the ideal cluster fractions according to Eq. ^ are used 
which arise from a purely stochastic distribution of the clus- 
ters. These jc^ do neither depend on the temperature nor the 
clusters' excess energies but are only determined by x and iij. 
The number W is then simply given by all possible arrange- 
ments of Np^ = xN A atoms and A^b = ( 1 ^ x)N B atoms on the 

= Na+Ns cation sites of the alloy, i.e.. 



^SRS _ 



A^aWb! 



(12) 



In the Stirling limit the mixing entropy in this case reduces to 

AS^^^ (x) = -ksN [x\nx + ( 1 - .^:) ln( 1 - x)] . (13) 

The ideal x'j according to Eq. ^ and A^^*^^ can be inter- 
preted as the high-temperature limit of the GQCA, since for 
increasing temperatures the xj^'^^ approach the x^ and, con- 
sequently, AS [cf. Eq. ©l approaches A^^'^^. 

(ii) The microscopic decomposition model (MDM). In this 
limiting case the cations of type A (B) are more likely to occur 
close to cations of type A (B). Consequently, only the clusters 
representing the two binary components AC and BC are al- 
lowed, with xM being the number of AC clusters and ( 1 —x)M 
the number of BC clusters. This is equivalent to merely a 
linear interpolation between the binary end components and, 
therefore, the results are equal to predictions from Vegard's 
rule^Sili. Furthermore, for positive excess energies Ae^ [cf. 
Eq. (|6]l] the x^*™ represent the low-temperature limit of the 
GQCA. The cluster fractions for the MDM are given by 



1 — X for j = 
X for j — J 
otherwise 



(14) 



The number of atomic configurations follows immediately 
from Eq. dTji by taking into account only the two clusters with 
j — and j — J, i.e.. 



W 



MDM 



Ml 



[xM]![(l-x)M]!" 



(15) 



Using the Xj we obtain for the mixing entropy 

A5'^™(x) = -/tBM[xlnx+ (1 -x)ln(l -x)] 

= -AS^^^ix), 
n 



(16) 



i.e., a reduced configurational entropy. In the MDM the result- 
ing length scale of decomposition of the alloy into the binary 
compounds is not well defined as usual in a thermodynamic 
approach. The minimum length scale is given by the assumed 
cluster size itself. Whether clusters of the same type (AC or 
BC) build larger regions of pure AC and BC material on a 
mesoscopic or macroscopic length scale depends on the clus- 
ter interaction (which is not taken into account). In general. 
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the MDM describes alloys that have been prepared under con- 
ditions where mixing does not lead to a gain of internal en- 
ergy. Indeed, in the MDM description one finds from Eq. ^ 
AU{x) = for the mixing contribution to the internal energy. 



III. QUANTUM MECHANICAL CALCULATIONS OF 
CLUSTERS 

A. Computational methods 

The statistical methods introduced in Section trace the 
properties of an alloy of composition x at temperature T back 
to the corresponding properties of individual clusters using 
Eq. (O. Therefore, the respective calculations have to be 
performed for each of the 22 (wz) and 16 (rs) clusters of 
Mg.,-Zni_ vO and Cd;i-Zni_ vO. 

Ground-state properties, such as the total energy of a cluster 
or its structural properties, are well described within density 
functional theory (DFT)'*^'*^ applying the supercell method, 
where a cluster is identified with the unit cell of an artificial 
crystal. Explicit calculations are performed using the Vienna 
Ab-initio Simulation Package (VASP)"*^^. For the semi-local 
generalized-gradient approximation (GGA) to the exchange 
and correlation (XC) potential in the Kohn-Sham equation^ 
we rely on the PW9 1 parametrization as suggested by Perdew 
and Wang'*^ ''^^. Although such a semi-local treatment of XC 
cannot correctly describe excitation properties (such as QP en- 
ergies or excitonic effects'*^'^^), it allows for a very good de- 
scription of lattice constants and cohesive energies. For the 
oxides studied in this work, the GGA tends to slightly un- 
derestimate binding energies and slightly overestimate lattice 
constants^^. However, these ground-state properties are not 
systematically improved by a different treatment of XC, for 
instance by adding an additional on-site Hubbard interaction 
U'^'^ or a certain amount of Fock exchange^". 

The electron-ion interaction is modeled within the 
projector-augmented wave (PAW) method^''^^. We treat the 
Mg2/7 and Mg3i as well as the 02i and 02p electrons as 
valence electrons. For Zn and Cd, the Zn 3d and Cd4d states 
have to be included as valence states, respectively, together 
with the Zn45 and CdSi electrons. The electronic wave func- 
tions are expanded into plane waves up to a cutoff energy of 
450 eV and the Brillouin zone (BZ) is sampled by 2x2x2 
(wz) or 4x4x4 (rs) Monkhorst-Pack (MP) k points^^. We 
checked that these parameters give results for the total ener- 
gies of the clusters that are converged within 0.08 %. This cor- 
responds to a very small level of inaccuracy, especially since 
the total energies of the different binary end components differ 
by at least 10 % of their total energy. 

In order to obtain the equilibrium values for the cell vol- 
umes, total energies, and bulk moduli we minimize the total 
energy for several cell volumes and fit the resulting volume 
dependence of the total energy to the Murnaghan equation of 
state^^. For each cluster geometry, we compute fully relaxed 
atomic positions, i.e., optimized lattice constants and internal 
cell parameters, ensuring that the forces acting on the ions are 
well below 5 meV/A. 



B. Results 

To describe the geometry of the two different crystal struc- 
tures on the same footing we discuss the average cation-cation 
(second-nearest neighbor) distance d2. The results for the 
cluster energies ej, the cation-cation distances d2j (averaged 
over all cation-cation distances of one cluster), the wz lattice 
parameters cj of the cluster cell and uj (averaged over all m's 
of one cluster), the cell volumes Vj, and the bulk moduli Bq j 
are listed in Table HIl for all wz cluster cells and Table Hill for 
all cluster classes of the rs crystal structure. 

The energies of the clusters Ej in Tables |ll] and [III] show 
an almost linear variation with the number rij of Mg (Cd) 
atoms, decreasing (increasing) from pure ZnO to pure MgO 
(CdO). Small negative (positive) non-linear deviations appear 
for Mgn^Zns-n^Og (Cd„^.Zn8_„^.08) with absolute values of a 
few meV (tenths of an eV). 

The cation-cation distances d2.j show a weak decrease 
(strong increase) with the number nj of Mg (Cd) atoms as 
expected from the corresponding values for the binary end 
components. While the nearest-neighbor distances or bond 
lengths (not listed in Tables HIl and HIB of the ZnO-MgO sys- 
tem decrease only very shghtly from dzn-o ~ 2.002/2.167 A 
to t/Mg-o ~ 1.997/2.128 A for wz/rs, we find a completely 
different situation for the ZnO-CdO system, where the bond 
lengths dcd-o ~ 2.230/2.388 A are rather different from 
dzn-o- For the bulk moduli Bo./ we observe a decrease with 
increasing number nj of Mg or Cd atoms, independent of the 
crystal structure. This trend follows the behavior of the co- 
valent radii of the cations that also show an anomalous trend 
of 1.36 A, 1.25 A, and 1.48 A along the row Mg, Zn, and 
Cd^^. For our calculated lattice parameters, equilibrium vol- 
umes, and bond lengths we find very good agreement with a 
similar study^. However, the deviations are remarkable for 
the total energies and, consequently, lead to different results 
for, e.g., the mixing free energy. We attribute these differ- 
ences to the use of the local-density approximation (LDA) to 
XC, and to the neglect of the Zn 3d and Cd 4-d electrons as 
valence electrons in Ref.[36l 



IV. THERMODYNAMIC PROPERTIES 

By including the temperature-dependent mixing entropy we 
also obtain a temperature dependence of the cluster fractions 
and, therefore, of the derived properties. In this work we 
discuss the alloy properties at two selected temperatures, (i) 
room temperature (T = 300 K) where many measurements 
are performed, and (ii) an exemplary growth temperature of 
r= HOOK. 



A. Excess energies 

The excess energy [cf. Eq. ©J Ae^ represents the energy 
of formation of a certain X„^Zn„_„^0„ cluster in the wz or rs 
structure with respect to ZnO and XO (X = Mg, Cd) in the re- 
spective crystal structure. They are fundamental quantities for 



6 




Cd fraction n./n Cd fraction n./n 



FIG. 2. (Color online) Excess energies per cation-anion pair Ae^ of the Mg„^Zn„_„^0„ (a), (c) and Cd„^Zn„_„jO„ (b), (d) clusters computed 
from Eq. The excess energies are shown for the rocksalt (red triangles) and the wurtzite (blue circles) crystal structures. In subfigures (a) 
and (b) [subfigures (c) and (d)] the Eq and £j of the same crystal structure [equilibrium crystal structure] have been used as level of reference. 
The lines linearly connect the end points vvz-ZnO and vvz-XO (solid blue) as well as r^-ZnO and rs-XO (dashed red), X = Mg,Cd. 



the GQCA, since they ultimately determine the cluster frac- 
tions [see Eq. ( fTTI )l and, therefore, the alloy statistics. In Fig.|2] 
the excess energies are plotted for both material combinations, 
i.e. ZnO-MgO and ZnO-CdO. Whereas in Figs. IS a) and (b) 
the ASj for wz and rs structure are shown, we visualize the 
energetic difference of these two lattice structures by plotting 
the ASj using the value of Eq for wz-ZnO and Sj for rs-XO in 
Figs. 12c) and (d). 

For MgZnO clusters [cf. Fig. |2|a)] we find negative excess 
energies (per anion-cation pair) with absolute values of less 
than 20 meV for both crystal structures. Contrary, for CdZnO 
clusters [cf. Fig.|2|b)] the excess energies are positive for both 
polymorphs and by one order of magnitude larger than the 
ones for MgZnO. As a consequence, in the isostructural limit 
with only a fourfold (wz) or a sixfold {rs) coordinated bonding 
configuration, we predict that mixing is possible for the ZnO- 
MgO system independent of the composition x, temperature 
T, and the boundary condition that fixes the cluster fractions. 
In the case of the ZnO-CdO system with its large difference of 
the bond lengths c/zn-o / "^Cd-o and, hence, a remarkable strain 
of a bond in the environment of the other cation species, the 
situation is different and will be discussed in detail in Sec- 
tion |IVC] 

Due to the changed reference energy of the excess ener- 
gies in Fig. |2jc) and their small absolute values in the case of 
ZnO-MgO, we find an almost linear variation between the end 



points. The crossing point near nj/n w 2/3 indicates that for 
small Mg fractions the fourfold coordination is more stable 
while for larger Mg fractions the sixfold one seems to be en- 
ergetically favored. This is in agreement with the occurrence 
of the cubic crystal structure above a composition of x « 0.67 
in PLD samples'**. However, there are also experimental— 
or other theoretical results '^'^^'^^ contradicting our findings. 
More precisely, Refs. and ,36i find for ri-ZnO and wz-ZnO 
almost the same energy difference as we do for rs-MgO and 
wz-MgO and vice versa. There are strong indications that the 
use of the LDA instead of the GGA can cause such devia- 
tions of the total energy differences. Our presumption is sup- 
ported by a comparison to other GGA calculations that con- 
firm our result of 0.3 eV as the energy difference between rs- 
ZnO and wz-ZnO as well as the 0. 15 eV between rs-MgO and 
wz-MgO^^^. LDA values that have been reported seem to 
lead to opposite differences for MgO and ZnO'^'^^'^^'^^'^^. In 
addition. Fan et alr'^ do not include the d electrons of Zn or 
Cd in their calculations which may further influence the rs-wz 
splittings. 

We observe a different situation for the ZnO-CdO system 
[cf. Fig. |2jd)], since the absolute values of the excess ener- 
gies are comparable to the energy differences of the rs and wz 
polymorphs. Consequently, the deviations from the linear in- 
terpolation between the end components are much larger. Fur- 
thermore, the wz and rs polymorphs of CdO are energetically 
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TABLE II. Cluster energies per cation-anion pair £j, cation-cation 
distances d2j, lattice parameters cj and uj, volume per cation-anion 
pair Vj, and bulk moduli Bq j for 16-atom clusters of Mg„^Zns_„^08 
(first line of each class) and Cd„^Zn8_„^08 (second line of each 
class) in wz structure. 



class j 


(eV/pair) 


dl.j 
(A) 


(A) 


Uj 


(A3) 


(GPa) 





-9.05 


3.270 


5.303 


0.378 


24.72 


141.0 


1 


-9.41 


3.268 


5.288 


0.380 


24.68 


137.5 




—8.87 


3.315 


5.368 


0.379 


25.75 


121.3 


2 


-9.76 


3.270 


5.275 


0.381 


24.73 


131.6 




-8.69 


3.359 


5.452 


0.378 


26.80 


110.2 


3 


-9.76 


3.266 


5.275 


0.381 


24.64 


131.3 




— 8.70 


3.357 


5.405 


0.382 


26.74 


112.7 


4 


-9.76 


3.266 


5.271 


0.381 


24.64 


131.1 




-8.71 


3.358 


5.445 


0.379 


26.77 


112.5 


5 


-10.12 


3.265 


5.243 


0.383 


24.61 


121.5 




-8.51 


3.408 


5.565 


0.374 


27.98 


114.2 


6 


-10.12 


3.265 


5.252 


0.382 


24.60 


121.6 




-8.55 


3.405 


5.508 


0.380 


27.91 


117.2 


7 


-10.12 


3.265 


5.260 


0.382 


24.60 


121.4 




—8.54 


3.404 


5.459 


0.384 


27.87 


116.0 


8 


-10.47 


3.264 


5.211 


0.385 


24.57 


118.4 




-8.33 


3.455 


5.705 


0.368 


29.16 


103.7 


9 


-10.47 


3.264 


5.235 


0.384 


24.58 


118.1 




-8.37 


3.451 


5.542 


0.383 


29.03 


104.0 


10 


-10.47 


3.264 


5.232 


0.384 


24.57 


118.0 




-8.38 


3.452 


5.591 


0.378 


29.07 


105.9 


11 


-10.47 


3.264 


5.246 


0.383 


24.58 


117.3 




-8.39 


3.447 


5.467 


0.390 


28.89 


104.1 


12 


-10.47 


3.263 


5.233 


0.384 


24.57 


118.8 




—8.41 


3.450 


5.583 


0.380 


29.03 


107.6 


13 


-10.47 


3.264 


5.241 


0.383 


24.57 


117.6 




-8.40 


3.449 


5.531 


0.384 


28.98 


105.1 


14 


-10.82 


3.261 


5.216 


0.385 


24.52 


116.8 




— 8.24 


3.495 


5.576 


0.387 


30.13 


95.9 


15 


-10.82 


3.261 


5.210 


0.385 


24.51 


116.7 




-8.25 


3.497 


5.636 


0.381 


30.22 


98.2 


16 


-10.82 


3.262 


5.198 


0.386 


24.53 


114.7 




-8.21 


3.501 


5.711 


0.374 


30.33 


96.3 


17 


-11.17 


3.260 


5.183 


0.387 


24.48 


113.2 




-8.11 


3.544 


5.714 


0.382 


31.48 


98.4 


18 


-11.17 


3.260 


5.191 


0.387 


24.49 


110.3 




-8.10 


3.542 


5.647 


0.387 


31.38 


96.6 


19 


-11.17 


3.261 


5.182 


0.387 


24.50 


109.8 




-8.08 


3.546 


5.727 


0.380 


31.51 


96.0 


20 


-11.52 


3.259 


5.158 


0.389 


24.46 


110.5 




-7.95 


3.592 


5.762 


0.384 


32.74 


92.4 


21 


-11.86 


3.258 


5.120 


0.391 


24.41 


104.4 




-7.82 


3.640 


5.841 


0.384 


34.10 


90.6 



much closer to each other than the corresponding ZnO poly- 
morphs. Merely plotting the excess energies [see Fig. |2d)] 
suggests a wide range up to rij/n 0.95 for the stability of 
wz which contradicts the majority of experiment a 1 2, 1 3^23,24,60 ^ 
Taking the cluster statistics into account we discuss this issue 
in more detail later [cf. Sections [IV Cl and llV Dl . 



TABLE III. Cluster energies per cation-anion pair e^, cation-cation 
distances d2,j, volume per cation-anion pair Vj, and bulk moduli Bq j 
for 16-atom clusters of Mg„^Zn8_„^,08 (first line of each class) and 
Cdn^ Zng-n^Og (second line of each class) in rs structure. 



class j 












(eV/pair) 


(A) 


(A3) 


(GPa) 





-8.75 


3.065 


20.36 


166.3 


1 


-9.16 


3.059 


20.23 


165.4 




-8.57 


3.111 


21.29 


162.1 


2 


-9.57 


3.052 


20.09 


162.7 




-8.42 


3.152 


22.14 


154.9 


3 


-9.57 


3.053 


20.11 


162.9 




-8.34 


3.165 


22.43 


155.4 


4 


-9.98 


3.044 


19.95 


160.9 




-8.31 


3.189 


22.92 


147.5 


5 


-9.98 


3.045 


19.97 


155.0 




-8.24 


3.202 


23.21 


147.3 


6 


-10.39 


3.037 


19.81 


158.5 




-8.22 


3.224 


23.69 


144.9 


7 


-10.39 


3.037 


19.80 


158.7 




-8.22 


3.224 


23.70 


144.9 


8 


-10.39 


3.038 


19.83 


159.0 




-8.15 


3.236 


23.95 


145.0 


9 


-10.39 


3.039 


19.85 


159.2 




-8.09 


3.247 


24.22 


144.3 


10 


-10.80 


3.030 


19.67 


156.7 




-8.09 


3.267 


24.66 


138.4 


1 1 
i 1 


iU. ou 


D.UJ i 


1 Q 70 

1 7. / U 






-8.03 


3.280 


24.95 


132.5 


12 


-11.20 


3.023 


19.54 


156.3 




-7.99 


3.315 


25.75 


114.1 


13 


-11.21 


3.024 


19.56 


157.4 




-8.00 


3.344 


26.39 


101.8 


14 


-11.61 


3.016 


19.40 


153.5 




-7.90 


3.343 


26.42 


130.8 


15 


-12.01 


3.009 


19.26 


152.0 




-7.84 


3.378 


27.25 


126.9 



B. Internal energy and entropy 

For a general discussion of the thermodynamic trends re- 
sulting from the three different statistics we compare the con- 
tributions AU [Eq. ©I and AS [Eq. ( fTOl iI to the mixing free 
energy AF as calculated from the GQCA with results obtained 
using the ideal weights, and the weights from the MDM for 
wz-MgTZni_,.0. 

From the mixing contribution to the internal energy 
AU{x, T) of wz-MgxZni_.vO it can be seen [cf. Fig. [3 a)] that 
for low temperatures (T = 30 K) only the clusters with the 
lowest energies contribute in the GQCA. For higher temper- 
atures (T = 300 K) the weights of clusters with higher ex- 
cess energies are larger. The SRS model finally represents the 
high-temperature limit of the GQCA in which all wz clusters 
contribute. The MDM does not lower the internal energy of 
the system. 

In Fig. [3jb) the corresponding results for the mixing en- 
tropy AS{x, T) are shown. Again, the MDM is the most triv- 
ial case with the smallest entropy term, since only two clus- 
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FIG. 3. (Color online) Mixing contribution to the internal energy 
MJ (x, T) (a) and mixing entropy A5(x, T) (b) for MgvZni_j-0 in wz 
structure. The solid curves are obtained using the GQCA cluster 
fractions for T = 30 K (blue) and T = 300 K (red). The dotted curves 
result from the MDM and the dash-dotted curves are obtained using 
the ideal cluster fractions. For comparison the blue dots represent the 
excess energies computed from Eq. l|6). All quantities are normalized 
to cation-anion pairs. 



ters contribute. The SRS model is obviously maximizing 
the entropy which is consistent with what we expect from a 
high-temperature limit, i.e., the high-temperature curve of the 
GQCA approaches the SRS limit. The low-temperature curve 
of the GQCA shows an interesting behavior, since the entropy 
in this case clearly features pronounced minima for composi- 
tions X where only single clusters substantially contribute. 



C. Mixing free energy 
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FIG. 4. (Color online) Mixing free energy AF(.x, T) of Mg^^Zni- 
(a) and CdvZni_j:0 (b) alloys with fixed sixfold {rs) or fourfold (wz) 
bonding configuration versus composition x. The solid curves are 
obtained using cluster fractions from the GQCA. The dotted curves 
result from the MDM and the dash-dotted curves are obtained using 
the ideal cluster fractions. In each case we depict results for T = 
300 K (blue curves) and T = 1 100 K (red curves). 



The Gibbs free energy is the thermodynamic potential 
which describes the equilibrium of a system for fixed tempera- 
ture and pressure. At low pressures of about 1 atmosphere the 
difference between the Gibbs free energy and the Helmholtz 
free energy vanishes. Therefore, the Helmholtz mixing free 
energy is used as the central thermodynamic quantity that gov- 
erns, at least for solids and low pressures, the equilibrium al- 
loy state and, therefore, the cluster fractions x^'^'^^. 

In Fig. |4]the mixing free energy AF{x^ T) for Mg^Zni_^0 
and Cd^Zni_rO alloys in the wz and the rs crystal structures 
is plotted versus compositions. For MgtZni_vO it is obvious 
that AF < for all x and T, which is, according to Eqs. (|4]i 
and (|5]l, a direct consequence of the negative excess energies 
for MgnjZnn-iijOn clusters [cf. Fig.|2a)]. For that reason all 
three statistics agree in predicting wz- and r5-Mg^Zni_iO to 
be a random alloy without a miscibility gap^^, independent 
of X and T. There is no tendency for binodal or spinodal de- 
composition of the alloy. The results using the GQCA cluster 
weights and the ones from the SRS model are almost identical, 
while the energy gain is smaller within the MDM. 

For Cd^Zn i the behavior of AF is qualitatively different 
as can be seen from Fig.|4jb). It shows that for T — 300 K and 



both crystal structures the GQCA results agree very well with 
the MDM, i.e., we find the alloy being almost entirely decom- 
posed into the binary clusters. For the higher temperature of 
r = 1100 K the GQCA result approaches the SRS curve for 
the wurtzitic alloy, especially for jc < 0.1 andx > 0.8. How- 
ever, for ri-Cd.vZni_.^0 at T = 1100 K the GQCA curve still 
agrees with the MDM. In this case the system has a strong ten- 
dency for decomposition into the binary materials and the SRS 
limit can only be reached at very high temperatures. Besides, 
for Cd^Zni_vO the AF curves versus x are rather asymmetric 
for the SRS and show minima and inflection points with po- 
sitions that strongly depend on the temperature and the local 
crystal structure. This is an indication of a phase transition 
between random and phase-separated alloys and will be dis- 
cussed as such later 

In the following, we compare to free energy curves that 
are computed using the energies of the respective equilibrium 
crystal structures (r^s-MgO, wz-ZnO, ri-CdO) as level of ref- 
erence, in contrast to the results in Fig.|4]where the energies of 
the binary components in the same crystal structure have been 
used. For Mg^Zni_^0 the mixing free energy curves for the 
wz alloys and the rs alloys in Fig.|3a) intersect atx w 0.67, in- 
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FIG. 5. (Color online) Mixing free energy AF(x,T) of Mg^-Zni_vO 
(a) and CdvZni_,0 (b) alloys versus composition x. The curves are 
obtained using the GQCA with wurtzite clusters only (solid), rocksalt 
clusters only (dash-dotted), or both types of clusters (dashed). In all 
cases results are shown for T = 300 K (blue curves) and T = 1 100 K 
(red curves). The equilibrium crystal structures wz (x = 0) and rs 
(x= I) have been used as energy zero. 



FIG. 6. (Color online) Mixing free energy AF{x,T) of MgtZni_ vO 
(a) and CdiZni_ vO (b) alloys versus composition x foi T = 300 K 
(blue) and T = 1100 K (red). The solid curves are computed using 
cluster fractions from the GQCA. The dotted curves are obtained for 
the MDM while the dash-dotted curves are calculated using the ideal 
cluster fractions. All curves result from the combined statistics with 
both wurtzite- and rocksalt-type clusters. The respective composition 
end points have been used as zero (see text). 



dependent of the temperature. We interpret this as a tendency 
for a transition from preferred fourfold coordination (wz) to 
preferred sixfold coordination (rs) at that composition under 
equihbrium conditions. In addition, the difference of the mix- 
ing free energy per cluster of the heterostructural alloys and 
that of the respective isostructural cases exceeds 25 meV, i.e. 

at room temperature, for approximately 0.10 < x < 0.98 
(0.28 <x< 0.93) at r = 300 K (T = 1100 K). The addi- 
tional degree of freedom of varying also the atomic coordi- 
nation leads to the lowest mixing free energy [see Fig.|5ja)]. 
Hence, for these values of x both, rs as well as wz clusters, 
significantly contribute to the alloy material. 

The same tendencies have been observed in experimental 
studies. Minemoto et alM- (Vashaei et al:^) found predom- 
inantly wz structure below x < 0.46 (x < 0.34) and mainly 
rs structure for x > 0.62 (0.65 < jc < 0.97). Ohtomo et alP 
report on the occurrence of an impurity phase above Mg con- 
centrations of X « 0.33 which roughly coincides with the point 
in Fig. |5la) where the high-temperature curve for wz struc- 
ture starts to deviate from the curve for the mixed statistics. 
Also the X-ray diffraction measurements of thin-film sam- 
ples by Bundesmann et aO^ reveal hexagonal wz structure for 
X < 0.53 and cubic rs structure for x > 0.67. Films grown 
by reactive electron beam evaporation^' lead to hexagonal 
MgtZni_;rO up to X = 0.51 and to cubic MgiZni_xO above 
X = 0.55. In general, the amount of fourfold- or sixfold- 
coordinated atoms cannot be directly derived from the mix- 
ing free energy and will be discussed in more detail in Sec- 
tion |IVD] 

In the case of the CdiZni_xO alloys, we study the mixing 
free energy in Fig. |5lb) where we observe a crossing of the 
curves for the isostructural alloys roughly at a Cd content of 
X « 0.95. Moreover, we find that the result from the mixed 
statistics differs less than 25 meV (per cluster) from the curve 
for the pure wz structure up to compositions x of about 0.17 



(0.59) for r = 300 K (r = 1 100 K). The reason for that be- 
havior is the small energy difference between the rs-CdO and 
the wz-CdO phase (see Tables HIl and HIIli. 

Experimental studies of CdxZni_vO show an ambivalent 
picture: While two groups report very low thermodynamic 
solubility limits^^ of only x ss 0.07 or phase separation at 
even lower Cd concentrations^'*, another group observed Cd 
concentrations up to x = 0.32 in samples produced by means 
of highly non-equilibrium MBE'^. Unfortunately, they have 
not tried for higher concentrations. More importantly, the 
wz crystal structure has also been found for plasma-enhanced 
MOCVD layers'^ up to x = 0.697. Their result is confirmed 
by Ishihara et alS''^ who report a transition from wz to rs struc- 
ture at X = 0.7 for films deposited by MOCVD. Apparently, 
the change of the crystal structure occurs at lower Cd con- 
centrations than we predicted from the intersection of the 
curves for wz and rs in Fig.jSjb). On the other hand, the high- 
temperature curve from the mixed statistics shows significant 
deviations from the pure wz statistics above Cd contents of 
about 0.7 which may explain the experimental findings^". 

For a detailed discussion of the mixed statistics we present 
in Fig. |6]the results for the mixing free energy obtained from 
the statistics with wz and rs clusters. In this plot we use again 
the equilibrium crystal structures for each material as level 
of reference for the mixing free energies at x = and x = 1 . 
Since the ideal cluster fractions x^ according to Eq. (|9]l do not 

depend on the cluster energies, both the sum of the x^ for all 
wz clusters and also of all rs clusters give the same total weight 
of 0.5. Consequently, at x = and x = 1 clusters of both 
crystal structures contribute equally, albeit, we want to take 
this situation into account due to possible non-equilibrium 
growth conditions. Therefore, we set the mixing free ener- 
gies at X = and x = 1 to zero for each curve resulting from 
the SRS model in Fig. |6] since otherwise, AF(x = 0,r) > 
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FIG. 7. (Color online) Crystal-structure fractions x"'~ (dashed curves) 
and x" (dash-dotted curves) of MgvZni_iO (a) and Cd^Zni_iO (b) 
alloys versus composition x. We calculated the curves using the 
cluster fractions from the GQCA for T = 300 K (blue curves) and 
T = 1100 K (red curves). The solid green lines are obtained using 
the ideal cluster fractions. For comparison the black curves show the 
result from the MDM for both crystal structures. 



and AF{x = l,T) > hold for temperatures T > K due to 
the weights of the SRS model. 

First of all, Fig.|6]points out that the GQCA and the MDM 
coincide for both materials at low temperatures. Contrary, the 
SRS model and the GQCA result in a remarkable difference 
independent of the temperature. This is not surprising, since 
the SRS model neglects the large energetic differences be- 
tween the two crystal structures by merely assigning the ideal 
weights to the clusters. While we find for Mg^Zni-^O, which 
shows negative excess energies, that the ideal cluster fractions 
are energetically favored over the ones resulting from GQCA 
or the MDM [cf. Fig.lMa)], the opposite is true for Cd;tZni_,,0 
[cf. Fig.|6tb)], where the excess energies are positive. 



D. Phase transitions and phase diagram 

For a better understanding of the structural composition of 
the alloys as a function of x and T we calculate the contribu- 
tions of wz and rs clusters to the mixed statistics results. For 
that purpose we define the wz character x^'^ and the rs charac- 
ter x" of the alloy as the sum over the corresponding weights, 
i.e.. 



21 



37 



Y^xj andx" = ^ xj, 



(17) 



=22 



with + x" = 1 . For given composition x, temperature 
T, and statistical model x"" (x") measures the relative con- 
tribution of clusters with fourfold (sixfold) atomic coordi- 
nation to the studied (non-decomposed) XxZni_^0 sample 
(X = Mg,Cd). The results for x""' and x" as a function of x 
and T are shown for both pseudobinary materials in Fig. |7] 
This figure clarifies how the dominating crystal structure in 
the alloy depends on the thermodynamic conditions as well as 
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FIG. 8. (Color online) Phase diagram computed for CdvZni_^0 us- 
ing the SRS model for the cluster statistics. We show binodal (black) 
and spinodal (blue) curves for wz clusters only (dashed curves), rs 
clusters only (dash-dotted curves), as well as for the mixed statis- 
tics including both sixfold and fourfold atomic coordination (solid 
curves). 



the respective cluster statistics (GQCA, SRS, MDM). As dis- 
cussed before, the SRS model for the cluster fractions gives 
rise to equal contributions of rs and wz clusters (see Fig. |7]). 
We confirm our discussion of the results for the mixing free 
energy by finding the GQCA curves close to the MDM re- 
sults for both materials at room temperature. As expected, 
higher preparation temperatures tend to move the intersection 
x"" = x" to larger Mg or Cd molar fractions x. More specif- 
ically, we find that point at about x = 0.5 [T = 300 K) and 
X = 0.6 (r = 1100 K) for Mg.vZni_.,0, while the tempera- 
ture dependence is more pronounced for CdvZni_vO where 
the intersection for T = 1100 K occurs at about x « 0.87. 
Consequently, the local crystal structure and bonding config- 
uration of the CdTZni-iO alloy depends much more on the 
actual growing conditions which explains the ambivalent ex- 
perimental findings for that material system i^^i'^i^^'^^i^" 

As discussed before, the occurrence of two pronounced 
minima as well as inflection points in the curve for the mixing 
free energy of CdvZni_,0 obtained from the SRS model (cf. 
Fig- 111) is an indication for a possible phase separation. Due 
to its negative excess energies [cf. Fig.|2ja)] the MgvZni_^0 
alloy does not show such a behavior. To study the conse- 
quences we construct the common tangent line to the free 
energy curves obtained from the SRS model for Cd;fZni_;fO 
alloys at several T. The resulting tangent points xi = xi{T) 
and X2 — X2{T) describe the boundaries of the miscibility gap 
and their variation with the temperature defines the binodal 
line in the T -x phase diagram shown in Fig. [H] From the 
maximum of the resulting T -x curves we find the critical pa- 
rameters Tcrit — 1030 K and Xcrit — 0.34 for the wz statistics as 
well as Tcrit ~ 1940 K and Xcrit — 0.4 for the mixed statistics. 
For the pure rs statistics the critical temperature is larger than 
2500 K. Inside the miscibility gap the alloy consists of a mix- 
ture of Cdx,Zni_;cjO and Cd;i2Zni_;c20 with the mixing free 
energy 



X2-X 



AF{xuT)- 



X2 -XI 

which is lower than AF{x, T). 



x — xi 

X2-X1 



AF(x2,r), (18) 
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FIG. 9. (Color online) Negative cohesive energy E^oh (in eV per 
pair) of MgvZni- tO (a) and Cdj:Zni_vO (b) alloys versus composi- 
tion X. The curves are obtained using only wurtzite clusters (dashed 
line), only rocksalt clusters (dash-dotted line), or both type of clusters 
(solid line). We use cluster fractions from the GQCA for T = 300 K 
(blue curves) and T = 1100 K (red curves), as well as ideal cluster 
fractions (green curves). The MDM result coincides with the one for 
from GQCA at room temperature. The cohesive energies are also 
shown for the clusters in rocksalt (red triangles) and wurtzite (blue 
circles) crystal structure. 
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FIG. 10. (Color online) Lattice parameter c (in A) of the wurtzite 
structure for Mgj:Zni_tO (a) and CdfZni_tO (b) alloys with clus- 
ters in wz crystal structure only versus composition x. The black 
curves represent the MDM and the green curves the ideal cluster frac- 
tions. The GQCA results are shown for T = 300 K (blue curve) and 
r = 1 100 K (red curve). For Mg.jZni_ vO the GQCA curves and the 
SRS results coincide. For Cdj^Zni_iO the low-temperature GQCA 
coincides with the MDM and the high-temperature GQCA with the 
result using ideal cluster fractions. Blue dots give the results for the 
individual clusters. 



Furthermore, the two inflection points x'^ = x\ (T) and ~ 
X2(r) of the AF{x,T) curve from the SRS model for 
CdtZni_^0 [cf. Fig. lUb)] define the spinodal curve in the 
phase diagram in Fig. [8] Inside the intervals xi < x < x[ and 
x'2 < X < X2 the alloy is metastable against local decomposi- 
tion due to an energy barrier Contrary, spontaneous decom- 
position into two alloys with compositions xi andx2 happens 
within the interval x\ < x < Xj.^^ The actual distribution of 
the clusters depends via xi (T) and X2(T) on the composition 
X and temperature T (cf. Fig. [8). 

The phase diagram in Fig.|8]suggests the speculation that the 
low solubilities of ri-CdO in wz-ZnO ^^'^^ can be explained 
within the SRS model. For alloy preparation conditions that 
lead to a stochastic occupation of the cation lattice sites. 
Fig. |8]confirms a large miscibility gap starting at x w 0.08 at 
r w 830 K (wz clusters only), T w 1450 K (rs and wz clusters), 
or r w 2070 K (rs clusters only). 

E. Cohesive energies 

From the cluster energies ej as given in Tables HI] and |lll] 
we compute the cohesive energies (per cation-anion pair) by 
subtracting the total energies of the spin-polarized atoms, gMg, 
£zn, £cd, and eo, from the ground state energy of the jih 
cluster, i.e., 

Ecoh.j = -£j + -^£x H + eo, (19) 

8 8 

with X = Mg, Cd. For the binary end components this leads to 
cohesive energies of £^011 = 10. 18 eV (ri-MgO), 7.21 eV (wz- 
ZnO), and 6.00 eV (r.s-CdO) which are close to experimental 
values of 10.26 eV, 7.52 eV, and 6.40 eV (see compilation in 



Ref.[ll. InFig.|9]the averaged cohesive energies [cf. Eq. ©I 
per cation-anion pair, as computed using the cluster weights 
from the different statistics, are plotted for Mg^Zni-vO and 
CdrZni_^0 alloys versus composition x. From Fig. |9|a) it 
immediately becomes clear that for MgjcZni_;fO the variation 
of the cohesive energy with the alloy composition x is almost 
linear. This is due to the fact that the variation between the dif- 
ferent crystal structures as well as the variation of the excess 
energies with x for fixed crystal structure (cf. Fig.|2]i both are 
small compared to the large difference of E^-oh between ZnO 
and MgO. For that reason also the temperature variation of the 
GQCA results is not visible. Overall, neither for the rs nor the 
wz crystal structure there is a visible dependence on the actual 
statistics that has been used. 

For CdvZni_;tO we find from Fig. |9lb) that the results 
within the GQCA and the SRS model differ. The curves from 
the MDM cannot be distinguished from the low-temperature 
GQCA results. For the wz crystal structure deviations result- 
ing from the different statistics are small but visible. As dis- 
cussed before, the mixed statistics gives rise to results that are 
close to those obtained when taking only wz clusters into ac- 
count. 



V. STRUCTURAL AND ELASTIC PROPERTIES 

A. Fourfold bonding coordination 

In Section HV] we found that especially for the Cd^Zni-^O 
system a significant fraction of the alloy locally shows the wz 
crystal structure over a large composition range. For that rea- 
son we investigate the behavior of the c lattice constant and the 
M parameter as characteristic parameters of the wz structure. 
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FIG. 11. (Color online) Internal cell parameter ii of the wurtzite 
structure for Mgx Zn i _ (a) and Cdj^Zn i _ (b) alloys derived from 
clusters in wz crystal structure only versus composition x. The black 
curves represent the MDM and the green curves the SRS model. The 
GQCA results are shown for T = 300 K (blue curve) and T = 1 100 K 
(red curve). For Mg,Zni_vO the GQCA curves and the ideal cluster 
fractions coincide. For CdvZni_vO the low-temperature GQCA co- 
incides with the MDM. Blue dots give the results for the individual 
clusters involved. 

Figure [Tot a) shows that in the case of Mgj:Zni_;cO the re- 
sults from the GQCA for both temperatures and the curve ob- 
tained using the ideal cluster fractions almost coincide. The 
two statistics give rise to slightly larger values for c than the 
MDM and, therefore, also Vegard's law, which we find to be 
violated. Around x = 0.5 the deviation may reach values of 
about 0.5 %. As can be seen in Fig. [TOt a) this leads to large 
errors when the composition of an alloy is determined by mea- 
suring the c-lattice constant and applying Vegard's rule. 

Contrary, for Cdj:Zni_j;0 [cf. Fig. [TOtb)l the curves com- 
puted with the GQCA cluster fractions for high and low 
temperatures represent the two limiting cases. While the 
low-temperature curve coincides with the MDM, the high- 
temperature curve is close to the results from the SRS model 
and the deviations between the two cases are small. 

A comparison of the values cj for the individual clusters of 
both material systems shows that the variation is much larger 
for Cd„^Zn„_„^0„ than for Mg„^Zn„_„^0„. Overall, the devi- 
ations from the MDM (which is equivalent to Vegard's law) 
are larger for Mg;fZni_;fO and we find a remarkable bowing 
of the curves which renders a purely linear interpolation ques- 
tionable. Contrary, the bowing is smaller for CdTZni_;fO de- 
spite the stronger deviations for individual clusters, e.g., with 
rij — 3, 4, or 5 Cd atoms. 

In Fig.[TT]we also study the influence of the cluster statis- 
tics on the internal cell parameter u of the hexagonal lattice, 
which describes the translation of the cation sublattice with 
respect to the anion sublattice along the c direction. Again, 
for MgvZni_vO there is a noticeable deviation of all statis- 
tics (that are almost equal) from the MDM result, with a re- 
markable bowing of the curves. The u values of the individ- 
ual clusters show a larger scattering for Cd„^Zn„_„^0„ than 
for Mg„^Zn„_„^0„. In the case of Cd;cZni_vO again the low- 
temperature GQCA curve coincides with the MDM result. 



FIG. 12. (Color online) Average cation-cation (second-nearest 
neighbor) distance d2 of Mg,Zni_iO (a) and Cd,vZni_,0 (b) al- 
loys versus composition x. The curves are obtained using the GQCA 
with wurtzite clusters only (dashed curves), rocksalt clusters only 
(dash-dotted curves), or both types of clusters (solid curves). For 
the GQCA we present results for T = 300 K (blue curves) and 
T = 1100 K (red curves). For comparison we show results ob- 
tained using the MDM (black curves) and the SRS model (green 
curves). The cation-cation distances are also depicted for the indi- 
vidual clusters in rocksalt (red triangles) and wurtzite (blue circles) 
crystal structure. 

while the high-temperature GQCA curve is close to the one 
obtained using ideal cluster fractions. 

B. Cation-cation distance 

Both crystal structures that are taken into account to de- 
scribe the pseudobinary alloy system, rs and wz, are char- 
acterized by different lattice parameters. While the cubic rs 
structure is described by only one lattice constant ao, there are 
three independent parameters a, c, and u for the hexagonal wz 
structure which renders a universal description of the lattice 
structure of the alloys more difficult. As mentioned before, we 
solve this problem by restricting our discussion to the average 
cation-cation distance d2{x). This quantity can be defined in 
a unique manner even for different crystal structures and is 
experimentally accessible at the same time. For the rs lattice 
d2 is simply related to the lattice constant qq by d2 — flo/\/2 
and for the ideal wz lattice with c/a = \/8/3 it holds d2 — a. 
We average over all cation-cation distances in one cluster to 
obtain the d2.j (cf. Tables and Ullb . Using the respective 
weights from the GQCA, the SRS model, and the MDM, the 
configurational average is calculated from Eq. (O and the re- 
sults are plotted in Fig.[T2]for Mgj:Zni_vO and CdxZni_;tO. 

For Mg;tZni_iO [cf. Fig.[T2l'a)l we find for the pure wz and 
pure rs cluster geometries that the three alloy statistics give 
almost the same results and, therefore, also agree well with 
Vegard's rule. For each fixed crystal structure the variation 
of d2 versus x is small, whereas d2 differs up to 8 % between 
the wz and the rs structure of a material. While the GQCA 
results from the mixed statistics agrees well with the MDM 
for T — 300 K, we find significant deviations for T = 1 100 K. 
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FIG. 13. (Color online) Bulk modulus Bo (in GPa) of Mg.iZni_, 0(a) 
and CdiZni_j:0 (b) alloys versus composition x. The curves are ob- 
tained using the GQCA with wurtzite clusters only (dashed curves), 
rocksalt clusters only (dash-dotted curves), or both type of clusters 
(solid curves). For the GQCA we present results for T = 300 K (blue 
curves) and T = 1 100 K (red curves). For comparison we also show 
results obtained using the MDM (black curves) and the SRS model 
(green curves). The bulk moduli are also depicted for the individual 
clusters in rocksalt (red triangles) and wurtzite (blue circles) crystal 
structure. 



In the high-temperature case the cation-cation distance in the 
pseudobinary alloy remains very close to the value of the 
pure wz structure up to compositions of about x w 0.2 [cf. 
Fig- El a)]- For higher Mg contents c/2 approaches the value 
of rs-MgO. 

For Cd;fZni_^0 [cf. Fig.[T2|b)] the variation of d2 is much 
larger than for MgTZni_vO since the cation-cation distance 
varies much more between ZnO and CdO. The results from 
all three statistics for the pure wz clusters and the pure rs clus- 
ters agree well with Vegard's rule except of a small deviation 
for the SRS model for ri-CdO. Again the most remarkable 
change is observed for the mixed cluster expansion. At room 
temperature the curve connects the wz and rs end components 
almost linearly. For the higher temperature of T = 1 100 K the 
cation-cation distance closely follows the wz trend up to Cd 
concentrations of more than 0.5. Overall, the different struc- 
tural configurations can be clearly distinguished at fixed x via 
the average cation-cation distance. When both lattice struc- 
tures can occur in the alloy we even find a more pronounced 
temperature dependence. 



C. Bulk modulus 

Besides the cation-cation distance we also study the bulk 
modulus Bo for Mg;iZni_;fO and Cdj:Zni_jcO using the three 
statistics evolving from the cluster weights of the GQCA, the 
SRS model, and the MDM. The results are plotted together 
with the values for the individual clusters in Fig. [13] 

For MgtZni_;fO [cf. Fig.[T3la)l in the pure rs or wz struc- 
ture the GQCA results show a temperature dependence and 
also differ from the linear interpolation of the MDM, whereas 
the SRS curve is very close to the high-temperature result of 



the GQCA. The reason for this behavior is the stronger devi- 
ation of the bulk moduli for the individual clusters from the 
linear interpolation. As can be seen in Fig. [T3la) for the rs 
structure we find a certain cluster with an especially low bulk 
modulus at X = 0.625 which, therefore, coincides with the 
minimum of the GQCA curves. For the GQCA results using 
the mixed statistics we find a strong temperature dependence, 
since the room-temperature curve closely follows the MDM 
curve, while the high-temperature result matches the pure wz 
curve up to jc w 0.2 and then approaches the pure rs curve at 
x=l. 

In the case of Cd,Zni_,0 [see Fig. flJlb)] we find both 
GQCA curves for the pure rs lattice structure close to the 
MDM result. This is only true for the low-temperature re- 
sult of the pure wz structure, whereas the high-temperature 
curve differs remarkably in this case. For both rs and wz struc- 
ture, the SRS model shows the largest deviation from the lin- 
ear interpolation. Again this can be attributed to the stronger 
variation of the values for Bo of the individual clusters with 
respect to the MDM curve with the strongest deviation occur- 
ring for rs at X — 0.75 [cf. Fig.[T3jb)]. For the mixed statis- 
tics for Cd;fZni_vO we find again that the low-temperature 
curve closely follows the linear interpolation, while the high- 
temperature curve shows pronounced deviations. At low con- 
centrations it follows the high-temperature GQCA curve of 
the pure wz structure and does not even fully approach the rs 
limit due to the very small energetical difference of these two 
polymorphs for CdO (cf. Fig.|7]). 



VI. SUMMARY AND CONCLUSIONS 

In summary, we studied pseudobinary A;tBi-.iC alloys 
whose binary end components crystallize in different crystal 
structures under ambient conditions. To study Mg^Zni_^0 
and Cd.vZni_^0 with rs (MgO, CdO) and wz (ZnO) equiUb- 
rium lattice structures we adopted the cluster expansion for 
wz and presented a corresponding expansion for 16-atom rs 
clusters. To include simultaneously locally different bonding 
configurations with either fourfold (wz) or sixfold (rs) coor- 
dination the GQCA has been extended correspondingly. Be- 
sides the GQCA which minimizes the mixing free energy with 
respect to the cluster fractions xj we took two other distribu- 
tions Xj into account to simulate also certain non-equilibrium 
growth procedures. We have combined these methods with 
ab initio density functional calculations using the semi-local 
GGA to treat exchange and correlation. This provides knowl- 
edge of structural and energetic properties for each of the clus- 
ters which we use together with the cluster fractions to per- 
form the configurational averages and, hence, to derive an av- 
eraged alloy property for given composition and temperature. 

Correspondingly, we find that the alloys' thermodynamics 
is determined by the local chemical bonding and the cluster 
statistics and, therefore, the conditions under which the alloy 
is prepared. The mixing free energies indicate a change from 
preferred wz crystal structure to preferred rs crystal struc- 
ture at about x « 0.675 for MgtZni_;i-0 and about x « 0.95 
for Cd;rZni_vO. We find that for Mg^Zni_;rO as well as 
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CdrZni_j^O the random alloy is always the most favorable. 
Only for Cd,Zni_,0 in the SRS limit we observe the occur- 
rence of a miscibility gap and binodal as well as spinodal de- 
composition ranges with the critical parameters Tait = 1030 K 
and Xcrit = 0.34 for pure wz structure, T^nt — 1940 K and 
■^crit = 0.4 for the mixed statistics, and a critical temperature 
above 2500 K for pure rs crystal structure. In any case we 
were able to calculate even the temperature- and composition- 
dependent amount of wz and rs clusters in the system. The 
resulting ambivalent picture which is strongly dependent on 
the actual conditions, is verified by a variety of experimental 
results. For both material systems there are experiments that 
agree with our predictions but also other experiments with de- 
viating findings, showing that also measured properties of al- 
loy films have to be related to growth conditions, deposition 
method, possible subsequent temper steps, etc. 

For structural properties such as cation-cation distances 
and bulk moduli we also find a significant dependence on 
the statistics that has been employed for the alloy modeling. 



In addition, the pronounced temperature dependence of the 
structural parameters and also the boundary conditions due to 
the growth mechanism should lead to distinct experimental 
observations of the alloy composition. 

ACKNOWLEDGMENTS 

We thank Lara K. Teles and R. Goldhahn for valuable 
discussions. The research leading to these results has re- 
ceived funding from the European Community's Seventh 
Framework Programme (FP7/2007-2013) under grant agree- 
ment No. 211956 and by the Deutsche Forschungsgemein- 
schaft (Project No. Be 1346/20-1). We acknowledge grants of 
computer time from the Supercomputerzentrum Stuttgart and 
the HPC-EUROPA2 project (project number: 228398) with 
the support of the European Commission -Capacities Area- 
Research Infrastructures. A. Schleife thanks the Carl-Zeiss- 
Stiftung for support. 



^ D. C. Reynolds, D. C. Look, and B. Jogai, Solid State Commun. 
99, 873 (1996). 

2 D. C. Reynolds, D. C. Look, B. Jogai, and T. C. Collins, Phys. 

Rev. B 56, 13753 (1997). 
^ M. Joseph, H. Tabata, and T. Kawai, Jpn. J. Appl. Phys. Part 2 38, 

L1205 (1999). 

tj. Ozgiir, Y. I. Alivov, C. Liu, A. Teke, M. A. Reshchikov, 

S. Dogan, V. Avrutin, S.-J. Cho, and H. Morkoc, J. Appl. Phys. 

98, 041301 (2004). 
5 D. S. Ginley and C. Bright, MRS Bull. 25, 58 (2000). 

S. Tsukazaki, A. Ohtomo, T. Kita, Y. Ohno, H. Ohno, and 

M. Kawasaki, Science 315, 1388 (2007). 
^ A. P Ramirez, Science 315, 1377 (2007). 

^ L. Schmidt-Mende and J. L. MacManus-Driscol, Materials Today 
10, 40 (2007). 

^ A. Ohtomo and A. Tsukazaki, Semicond. Sci. Technol. 20, SI 

(2005) . 

R. Schmidt, B. Rheinlander, M. Schubert, D. Spemann, T. Butz, 
J. Lenzner, E. M. Kaidashev, M. Lorenz, A. Rahm, H. C. Sem- 
melhack, and M. Grundmann, Appl. Phys. Lett. 82, 2260 (2003). 
" S. Sadofev, S. Blumstengel, J. Cui, J. Puis, S. Rogaschewski, 
P. Schafer, Y. G. Sadofyev, and F. Henneberger, Appl. Phys. Lett. 
87, 091903 (2005). 

S. Sadofev, S. Blumstengel, J. Cui, J. Puis, S. Rogaschewski, 
P Schafer, and F. Henneberger, Appl. Phys. Lett. 89, 201907 

(2006) . 

S. Shigemori, A. Nakamura, J. Ishihara, T. Aoli, and J. Temmyo, 
Jpn. J. Appl. Phys. 43, L1088 (2004). 

W. Martienssen and H. Warlimont, eds., "Springer handbook of 
condensed matter and materials data," (Springer- Verlag, Berlin, 
2005). 

A. Schleife, F. Fuchs, J. FurthmuUer, and F. Bechstedt, Phys. Rev 
B 73, 245212 (2006). 

M. Sanati, G. L. W. Hart, and A. Zunger, Phys. Rev. B 68, 155210 
(2003). 

1^ E. Segnit and A. Holland, J. Am. Ceram. Soc. 48, 409 (1965). 
C. Bundesmann, A. Rahm, M. Lorenz, M. Grundmann, and 
M. Schubert, J. Appl. Phys. 99, 113504 (2006). 



12 



14 



22 



Z. Liu, Z. Mei, T. Zhang, Y. Liu, Y. Guo, X. Du, A. Hallen, J. Zhu, 
and A. Kuznetsov, |Journare Crystal Growth] 311, 4356 (2009), 
ISSN 0022-0248. 
20 X. Du, Z. Mei, Z. Liu, Y. Guo, T. Zhang, Y. Hou, Z. Zhang, 
Q. Xue, and A. Y. Kuznetsov, Advanced Materials 21, 4625 
(2009). 

2' J. Chen, W. Z. Shen, N. B. Chen, D. J. Qiu, and H. Z. Wu, Journal 
of Physics: Condensed Matter 15, L475 (2003). 
N. B. Chen, H. Z. Wu, D. J. Qiu, T. N. Xu, J. Chen, and W. Z. 
Shen, J. Phys.: Condens. Matter 16, 2973 (2004). 
T. Makino, Y. Segawa, M. Kawasaki, A. Ohtomo, R. Shiroki, 
K. Tamura, T. Yasuda, and H. Koinuma, Appl. Phys. Lett. 78, 
1237 (2001). 

2^ F. Bertram, S. Giemsch, D. Forster, J. Christen, R. Kling, 
C. Kirchner, and A. Waag, Appl. Phys. Lett. 88, 061915 (2006). 
W. R. Lambrecht, S. Limpigumnong, and B. Segall, MRS Internet 
J. Nitride Semicond. Res. 451, G6.8 (1999). 

2^ R. Thangavel, M. Rajagopalan, and J. Kumar, Solid State Com- 
mun. 137, 507 (2006). 

27 Y.-G. Kim, E.-C. Lee, and K. J. Chang, J. Korean Phys. Soc. 39, 
S92 (2001). 

2^ S.-H. Wei, L. G. Ferreira, J. E. Bernard, and A. Zunger, 

Phys. Rev. B|42, 9622 (1990). 
2' A.-B. Chen and A. Sher, Semiconductor Alloys (Plenum, New 

York, 1995). 

^" L. K. Teles, J. Furthmiiller, L. M. R. Scolfaro, J. R. Leite, and 
F Bechstedt, Phys. Rev. B 62, 2475 (2000). 
J. M. Sanchez, F. Ducastelle, and D. Gratias, Physica A 128, 334 
(1984). 

■^2 A. Zunger, in Statistics and Dynamics of Alloy Phase Transfor- 
mations, edited by P. E. A. Turchi and A. Gonis (Plenum Press, 
New York, 1994) Chap. First-principles statistical mechanics of 
semiconductor alloys and intermetallic compounds, p. 361. 



J. Z. Liu and A. Zunger, Phys. Rev. B 77, 205201 (2008). 



33 

^'^ A. Sher, M. van Schilfgaarde, A.-B. Chen, and W. Chen, Phys. 

Rev. B 36, 4279 (1987). 
35 C. Caetano, L. K. Teles, M. Marques, A. D. P. Jr., and L. G. Fer- 
reira, Phys. Rev. 74, 045215 (2006). 



15 



3^ X. F. Fan, H. D. Sun, Z. X. Shen, J.-L. Kuo, and Y. M. Lu, 
Journal of Physics: Condensed Matter 20, 235221 (2008). 
1. V. Maznichenko, A. Ernst, M. Bouhassoune, J. Henk, M. Dane, 
M. Liiders, P. Bruno, W. Hergert, 1. Mertig, Z. Szotek, and W. M. 
Temmerman, Phys. Rev. B 80, 144101 (2009). 

3^ J. W. D. Connolly and A. R. Williams, Ph ys. Rev. B| 27, 5169 
(1983). 

A. Zunger, S.-H. Wei, L. G. Ferreira, and J. E. Bernard, 
Phys. Rev. Lett. 65, 353 (1990). 

40 S.-H. Jeong, B.-S. Kim, and B.-T. Lee, Applied Physics Letters] 
82, 2625 (2003). 

41 L. Vegard, Zeitschrift fur Physik| 5, 17 (1921). 

42 p Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

43 W. Kohn and L. J. Sham, Phys. Rev^ 140 , A1133 (1965). 

44 G. Kresse and J. Furthmuller, PhysTRev . B 54, 1 1 169 (10 1996). 
4^ J. P. Perdew, in Electronic Structure of Solids '91, edited by P. Zi- 

esche and H. Eschrig (Akademie-Verlag, Berlin, 1991) p. 11. 
46 J. P Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992). 
4^ A. Schleife, F. Fuchs, C. Rodl, J. Furthmuller, and F. Bechstedt, 
physica status solidi (b) 246, 2150 (2009). 

A. Schleife, C. Rodl, F. Fuchs, J. Furthmuller, and F. Bechstedt, 

|Phys. Rev. B| 80, 035112 (2009). 

49 A. Janotti, D. Segev, and C. G. Van de Walle, [PhyTRev. B| 74, 

045202 (2006). 

^" F. Oba, A. Togo, 1. Tanaka, J. Paier, and G. Kresse, [Phys. Rev."^ 
77, 245202 (2008). 



48 



51 P E. Blochl, Phys. Rev. B 50, 17953 (1994). 

52 G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

53 H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976). 

54 F D. Murnaghan, Proc. Natl. Acad, of Sci. USA 30, 244 (1944). 

55 Table of Periodic Properties of Elements (Sargent- Welch, Skokie 
Illinois, 1980). 

56 L. Zhang and H. Huang, [Applied Physics Letters] 90, 023115 
(2007). 

5^ J. Carrasco, F. lllas, and S. T. Bromley, [Phys. Rev. LeE] 99, 
235502 (2007). 

5** A. Seko, F. Oba, A. Kuwabara, and 1. Tanaka, iRiys. Rev. B[ 72, 
024107 (2005). 

59 S. Limpijumnong and S. Jungthawan, [Phys. Rev. B[ 70, 054104 

(2004) . 

6" J. Ishihara, A. Nakamura, S. Shigemori, T. Aoki, and J. Temmyo, 

Applied Physics Letters 89, 091914 (2006). 
6' T. Minemoto, T. Negami, S. Nishiwaki, H. Takakura, and 

Y. Hamakawa, Thin Solid Films 372, 173 (2000). 

62 Z. Vashaei, T. Minegish i, H. S uzuki, T. Hanada, M. W. Cho, 
T. Yao, and A. Setiawan, [Journal of Applied Physics| 98, 054911 

(2005) . 

63 A. Ohtomo, M. Kawasaki, T. Koida, K. Masubuchi, H. Koniuma, 
Y. Sakurai, Y. Yoshida, T. Yasuda, and Y. Segawa, Appl. Phys. 
Lett. 72, 2466 (1998). 



